Substantial viral diversity in bats and rodents from East Africa: insights into evolution, recombination, and cocirculation.

BACKGROUND: Zoonotic viruses cause substantial public health and socioeconomic problems worldwide. Understanding how viruses evolve and spread within and among wildlife species is a critical step when aiming for proactive identification of viral threats to prevent future pandemics. Despite the many proposed factors influencing viral diversity, the genomic diversity and structure of viral communities in East Africa are largely unknown. RESULTS: Using 38.3 Tb of metatranscriptomic data obtained via ultradeep sequencing, we screened vertebrate-associated viromes from 844 bats and 250 rodents from Kenya and Uganda collected from the wild. The 251 vertebrate-associated viral genomes of bats (212) and rodents (39) revealed the vast diversity, host-related variability, and high geographic specificity of viruses in East Africa. Among the surveyed viral families, Coronaviridae and Circoviridae showed low host specificity, high conservation of replication-associated proteins, high divergence among viral entry proteins, and frequent recombination. Despite major dispersal limitations, recurrent mutations, cocirculation, and occasional gene flow contribute to the high local diversity of viral genomes. CONCLUSIONS: The present study not only shows the landscape of bat and rodent viromes in this zoonotic hotspot but also reveals genomic signatures driven by the evolution and dispersal of the viral community, laying solid groundwork for future proactive surveillance of emerging zoonotic pathogens in wildlife. Video Abstract.

Genomic characterization of human respiratory syncytial virus circulating in Islamabad, Pakistan, during an outbreak in 2022-2023.

In this study, conducted at the National Institute of Health, Islamabad, during an outbreak of human respiratory syncytial virus (hRSV) from December 2022 to January 2023, the first whole-genome sequences of hRSV isolates from Islamabad, Pakistan, were determined. Out of 10 positive samples, five were sequenced, revealing the presence of two genotypes: RSV-A (GA2.3.5, ON1 strain) and RSV-B (GB5.0.5.a, BA-10 strain). A rare non-synonymous substitution (E232G) in G the protein and N276S in the F protein were found in RSV-A. In RSV-B, the unique mutations K191R, Q209R, and I206M were found in the F protein. These mutations could potentially influence vaccine efficacy and viral pathogenicity. This research underscores the importance of genomic surveillance for understanding RSV diversity and guiding public health responses in Pakistan.

A giant virus infecting the amoeboflagellate Naegleria.

Giant viruses (Nucleocytoviricota) are significant lethality agents of various eukaryotic hosts. Although metagenomics indicates their ubiquitous distribution, available giant virus isolates are restricted to a very small number of protist and algal hosts. Here we report on the first viral isolate that replicates in the amoeboflagellate Naegleria. This genus comprises the notorious human pathogen Naegleria fowleri, the causative agent of the rare but fatal primary amoebic meningoencephalitis. We have elucidated the structure and infection cycle of this giant virus, Catovirus naegleriensis (a.k.a. Naegleriavirus, NiV), and show its unique adaptations to its Naegleria host using fluorescence in situ hybridization, electron microscopy, genomics, and proteomics. Naegleriavirus is only the fourth isolate of the highly diverse subfamily Klosneuvirinae, and like its relatives the NiV genome contains a large number of translation genes, but lacks transfer RNAs (tRNAs). NiV has acquired genes from its Naegleria host, which code for heat shock proteins and apoptosis inhibiting factors, presumably for host interactions. Notably, NiV infection was lethal to all Naegleria species tested, including the human pathogen N. fowleri. This study expands our experimental framework for investigating giant viruses and may help to better understand the basic biology of the human pathogen N. fowleri.

Evaluation of different genomic regions of rotavirus B and rotavirus C for development of real-time RT‒PCR assays.

BACKGROUND: The causative agents of diarrhea, rotavirus B (RVB) and rotavirus C (RVC) are common in adults and patients of all age groups, respectively. Due to the Rotavirus A (RVA) vaccination program, a significant decrease in the number of gastroenteritis cases has been observed globally. The replacement of RVA infections with RVB, RVC, or other related serogroups is suspected due to the possibility of reducing natural selective constraints due to RVA infections. The data available on RVB and RVC incidence are scant due to the lack of cheap and rapid commercial diagnostic assays and the focus on RVA infections. The present study aimed to develop real-time RT‒PCR assays using the data from all genomic RNA segments of human RVB and RVC strains available in the Gene Bank. RESULTS: Among the 11 gene segments, NSP3 and NSP5 of RVB and the VP6 gene of RVC were found to be suitable for real-time RT‒PCR (qRT‒PCR) assays. Fecal specimens collected from diarrheal patients were tested simultaneously for the presence of RVB (n = 192) and RVC (n = 188) using the respective conventional RT‒PCR and newly developed qRT‒PCR assays. All RVB- and RVC-positive specimens were reactive in their respective qRT‒PCR assays and had Ct values ranging between 23.69 and 41.97 and 11.49 and 36.05, respectively. All known positive and negative specimens for other viral agents were nonreactive, and comparative analysis showed 100% concordance with conventional RT‒PCR assays. CONCLUSIONS: The suitability of the NSP5 gene of RVB and the VP6 gene of RVC was verified via qRT‒PCR assays, which showed 100% sensitivity and specificity. The rapid qRT‒PCR assays developed will be useful diagnostic tools, especially during diarrheal outbreaks for testing non-RVA rotaviral agents and reducing the unnecessary use of antibiotics.

Early detection of emerging viral variants through analysis of community structure of coordinated substitution networks.

The emergence of viral variants with altered phenotypes is a public health challenge underscoring the need for advanced evolutionary forecasting methods. Given extensive epistatic interactions within viral genomes and known viral evolutionary history, efficient genomic surveillance necessitates early detection of emerging viral haplotypes rather than commonly targeted single mutations. Haplotype inference, however, is a significantly more challenging problem precluding the use of traditional approaches. Here, using SARS-CoV-2 evolutionary dynamics as a case study, we show that emerging haplotypes with altered transmissibility can be linked to dense communities in coordinated substitution networks, which become discernible significantly earlier than the haplotypes become prevalent. From these insights, we develop a computational framework for inference of viral variants and validate it by successful early detection of known SARS-CoV-2 strains. Our methodology offers greater scalability than phylogenetic lineage tracing and can be applied to any rapidly evolving pathogen with adequate genomic surveillance data.

Unveiling hidden structural patterns in the SARS-CoV-2 genome: Computational insights and comparative analysis.

SARS-CoV-2, the causative agent of COVID-19, is known to exhibit secondary structures in its 5' and 3' untranslated regions, along with the frameshifting stimulatory element situated between ORF1a and 1b. To identify additional regions containing conserved structures, we utilized a multiple sequence alignment with related coronaviruses as a starting point. We applied a computational pipeline developed for identifying non-coding RNA elements. Our pipeline employed three different RNA structural prediction approaches. We identified forty genomic regions likely to harbor structures, with ten of them showing three-way consensus substructure predictions among our predictive utilities. We conducted intracomparisons of the predictive utilities within the pipeline and intercomparisons with four previously published SARS-CoV-2 structural datasets. While there was limited agreement on the precise structure, different approaches seemed to converge on regions likely to contain structures in the viral genome. By comparing and combining various computational approaches, we can predict regions most likely to form structures, as well as a probable structure or ensemble of structures. These predictions can be used to guide surveillance, prophylactic measures, or therapeutic efforts. Data and scripts employed in this study may be found at https://doi.org/10.5281/zenodo.8298680.

Genomic analyses of Symbiomonas scintillans show no evidence for endosymbiotic bacteria but does reveal the presence of giant viruses.

Symbiomonas scintillans Guillou et Chrétiennot-Dinet, 1999 is a tiny (1.4 µm) heterotrophic microbial eukaryote. The genus was named based on the presence of endosymbiotic bacteria in its endoplasmic reticulum, however, like most such endosymbionts neither the identity nor functional association with its host were known. We generated both amplification-free shotgun metagenomics and whole genome amplification sequencing data from S. scintillans strains RCC257 and RCC24, but were unable to detect any sequences from known lineages of endosymbiotic bacteria. The absence of endobacteria was further verified with FISH analyses. Instead, numerous contigs in assemblies from both RCC24 and RCC257 were closely related to prasinoviruses infecting the green algae Ostreococcus lucimarinus, Bathycoccus prasinos, and Micromonas pusilla (OlV, BpV, and MpV, respectively). Using the BpV genome as a reference, we assembled a near-complete 190 kbp draft genome encoding all hallmark prasinovirus genes, as well as two additional incomplete assemblies of closely related but distinct viruses from RCC257, and three similar draft viral genomes from RCC24, which we collectively call SsVs. A multi-gene tree showed the three SsV genome types branched within highly supported clades with each of BpV2, OlVs, and MpVs, respectively. Interestingly, transmission electron microscopy also revealed a 190 nm virus-like particle similar the morphology and size of the endosymbiont originally reported in S. scintillans. Overall, we conclude that S. scintillans currently does not harbour an endosymbiotic bacterium, but is associated with giant viruses.

Nanopore adaptive sampling of a metagenomic sample derived from a human monkeypox case.

In 2022, a series of human monkeypox cases in multiple countries led to the largest and most widespread outbreak outside the known endemic areas. Setup of proper genomic surveillance is of utmost importance to control such outbreaks. To this end, we performed Nanopore (PromethION P24) and Illumina (NextSeq. 2000) Whole Genome Sequencing (WGS) of a monkeypox sample. Adaptive sampling was applied for in silico depletion of the human host genome, allowing for the enrichment of low abundance viral DNA without a priori knowledge of sample composition. Nanopore sequencing allowed for high viral genome coverage, tracking of sample composition during sequencing, strain determination, and preliminary assessment of mutational pattern. In addition to that, only Nanopore data allowed us to resolve the entire monkeypox virus genome, with respect to two structural variants belonging to the genes OPG015 and OPG208. These SVs in important host range genes seem stable throughout the outbreak and are frequently misassembled and/or misannotated due to the prevalence of short read sequencing or short read first assembly. Ideally, standalone standard Illumina sequencing should not be used for Monkeypox WGS and de novo assembly, since it will obfuscate the structure of the genome, which has an impact on the quality and completeness of the genomes deposited in public databases and thus possibly on the ability to evaluate the complete genetic reason for the host range change of monkeypox in the current pandemic.

Viral genome sequencing to decipher in-hospital SARS-CoV-2 transmission events.

The SARS-CoV-2 pandemic has highlighted the need to better define in-hospital transmissions, a need that extends to all other common infectious diseases encountered in clinical settings. To evaluate how whole viral genome sequencing can contribute to deciphering nosocomial SARS-CoV-2 transmission 926 SARS-CoV-2 viral genomes from 622 staff members and patients were collected between February 2020 and January 2021 at a university hospital in Munich, Germany, and analysed along with the place of work, duration of hospital stay, and ward transfers. Bioinformatically defined transmission clusters inferred from viral genome sequencing were compared to those inferred from interview-based contact tracing. An additional dataset collected at the same time at another university hospital in the same city was used to account for multiple independent introductions. Clustering analysis of 619 viral genomes generated 19 clusters ranging from 3 to 31 individuals. Sequencing-based transmission clusters showed little overlap with those based on contact tracing data. The viral genomes were significantly more closely related to each other than comparable genomes collected simultaneously at other hospitals in the same city (n = 829), suggesting nosocomial transmission. Longitudinal sampling from individual patients suggested possible cross-infection events during the hospital stay in 19.2% of individuals (14 of 73 individuals). Clustering analysis of SARS-CoV-2 whole genome sequences can reveal cryptic transmission events missed by classical, interview-based contact tracing, helping to decipher in-hospital transmissions. These results, in line with other studies, advocate for viral genome sequencing as a pathogen transmission surveillance tool in hospitals.

Phylogenetic lineage dynamics of global parainfluenza virus type 3 post-COVID-19 pandemic.

During the coronavirus disease 2019 (COVID-19) pandemic, outbreaks of parainfluenza virus type 3 (PIV-3) decreased due to infection control measures. However, a post-pandemic resurgence of PIV-3 has recently been observed. Nonetheless, the role of viral genetic epidemiology, possibly influenced by a genetic bottleneck effect, remains unexplored. We investigated the phylogenetic structure of the publicly available PIV-3 whole-genome and hemagglutinin-neuraminidase (HN) gene sequences spanning the last 65 years, including the COVID-19 pandemic. Sequences were retrieved from the nucleotide database of the National Center for Biotechnology Information using the search term "Human respirovirus 3." Sequence subsets covering all six genes of PIV-3 or the HN gene were designated as the whole-genome and HN surveillance data sets, respectively. Using these data sets, we constructed maximum-likelihood phylogenetic trees and performed a time-scaled analysis using a Bayesian SkyGrid coalescent prior. A total of 455 whole-genome and 1,139 HN gene sequences were extracted, revealing 10 and 11 distinct lineages, respectively, with >98% concurrence in lineage assignments. During the 2020 COVID-19 pandemic, only three single-lineage clusters were identified in Japan, Korea, and the USA. The inferred year of origin for PIV-3 was 1938 (1903-1963) for the whole-genome data set and 1955 (1930-1963) for the HN gene data set. Our study suggests that PIV-3 epidemics in the post-COVID era are likely influenced by a pandemic-driven bottleneck phenomenon and supports previous hypotheses suggesting s that PIV-3 originated during the early half of the 20th century.IMPORTANCEUsing publicly available parainfluenza virus type 3 (PIV-3) whole-genome sequences, we estimated that PIV-3 originated during the 1930s, consistent with previous hypotheses. Lineage typing and time-scaled phylogenetic analysis revealed that PIV-3 experienced a bottleneck phenomenon in Korea and the USA during the coronavirus disease 2019 pandemic. We identified the conservative hemagglutinin-neuraminidase gene as a viable alternative marker in long-term epidemiological studies of PIV-3 when whole-genome analysis is limited.

Conserved structures and dynamics in 5'-proximal regions of Betacoronavirus RNA genomes.

Betacoronaviruses are a genus within the Coronaviridae family of RNA viruses. They are capable of infecting vertebrates and causing epidemics as well as global pandemics in humans. Mitigating the threat posed by Betacoronaviruses requires an understanding of their molecular diversity. The development of novel antivirals hinges on understanding the key regulatory elements within the viral RNA genomes, in particular the 5'-proximal region, which is pivotal for viral protein synthesis. Using a combination of cryo-electron microscopy, atomic force microscopy, chemical probing, and computational modeling, we determined the structures of 5'-proximal regions in RNA genomes of Betacoronaviruses from four subgenera: OC43-CoV, SARS-CoV-2, MERS-CoV, and Rousettus bat-CoV. We obtained cryo-electron microscopy maps and determined atomic-resolution models for the stem-loop-5 (SL5) region at the translation start site and found that despite low sequence similarity and variable length of the helical elements it exhibits a remarkable structural conservation. Atomic force microscopy imaging revealed a common domain organization and a dynamic arrangement of structural elements connected with flexible linkers across all four Betacoronavirus subgenera. Together, these results reveal common features of a critical regulatory region shared between different Betacoronavirus RNA genomes, which may allow targeting of these RNAs by broad-spectrum antiviral therapeutics.

Cryptic proteins translated from deletion-containing viral genomes dramatically expand the influenza virus proteome.

Productive infections by RNA viruses require faithful replication of the entire genome. Yet many RNA viruses also produce deletion-containing viral genomes (DelVGs), aberrant replication products with large internal deletions. DelVGs interfere with the replication of wild-type virus and their presence in patients is associated with better clinical outcomes. The DelVG RNA itself is hypothesized to confer this interfering activity. DelVGs antagonize replication by out-competing the full-length genome and triggering innate immune responses. Here, we identify an additionally inhibitory mechanism mediated by a new class of viral proteins encoded by DelVGs. We identified hundreds of cryptic viral proteins translated from DelVGs. These DelVG-encoded proteins (DPRs) include canonical viral proteins with large internal deletions, as well as proteins with novel C-termini translated from alternative reading frames. Many DPRs retain functional domains shared with their full-length counterparts, suggesting they may have activity during infection. Mechanistic studies of DPRs derived from the influenza virus protein PB2 showed that they poison replication of wild-type virus by acting as dominant-negative inhibitors of the viral polymerase. These findings reveal that DelVGs have a dual inhibitory mechanism, acting at both the RNA and protein level. They further show that DPRs have the potential to dramatically expand the functional proteomes of diverse RNA viruses.

Bipartite viral RNA genome heterodimerization influences genome packaging and virion thermostability.

Multi-segmented viruses often multimerize their genomic segments to ensure efficient and stoichiometric packaging of the correct genetic cargo. In the bipartite Nodaviridae family, genome heterodimerization is also observed and conserved among different species. However, the nucleotide composition and biological function for this heterodimer remain unclear. Using Flock House virus as a model system, we developed a next-generation sequencing approach ("XL-ClickSeq") to probe heterodimer site sequences. We identified an intermolecular base-pairing site which contributed to heterodimerization in both wild-type and defective virus particles. Mutagenic disruption of this heterodimer site exhibited significant deficiencies in genome packaging and encapsidation specificity to viral genomic RNAs. Furthermore, the disruption of this intermolecular interaction directly impacts the thermostability of the mature virions. These results demonstrate that the intermolecular RNA-RNA interactions within the encapsidated genome of an RNA virus have an important role on virus particle integrity and thus may impact its transmission to a new host.IMPORTANCEFlock House virus is a member of Nodaviridae family of viruses, which provides a well-studied model virus for non-enveloped RNA virus assembly, cell entry, and replication. The Flock House virus genome consists of two separate RNA molecules, which can form a heterodimer upon heating of virus particles. Although similar RNA dimerization is utilized by other viruses (such as retroviruses) as a packaging mechanism and is conserved among Nodaviruses, the role of heterodimerization in the Nodavirus replication cycle is unclear. In this research, we identified the RNA sequences contributing to Flock House virus genome heterodimerization and discovered that such RNA-RNA interaction plays an essential role in virus packaging efficiency and particle integrity. This provides significant insight into how the interaction of packaged viral RNA may have a broader impact on the structural and functional properties of virus particles.

Functional properties of measles virus proteins derived from a subacute sclerosing panencephalitis patient who received repeated remdesivir treatments.

Subacute sclerosing panencephalitis (SSPE) is a rare but fatal late neurological complication of measles, caused by persistent measles virus (MeV) infection of the central nervous system. There are no drugs approved for the treatment of SSPE. Here, we followed the clinical progression of a 5-year-old SSPE patient after treatment with the nucleoside analog remdesivir, conducted a post-mortem evaluation of the patient's brain, and characterized the MeV detected in the brain. The quality of life of the patient transiently improved after the first two courses of remdesivir, but a third course had no further clinical effect, and the patient eventually succumbed to his condition. Post-mortem evaluation of the brain displayed histopathological changes including loss of neurons and demyelination paired with abundant presence of MeV RNA-positive cells throughout the brain. Next-generation sequencing of RNA isolated from the brain revealed a complete MeV genome with mutations that are typically detected in SSPE, characterized by a hypermutated M gene. Additional mutations were detected in the polymerase (L) gene, which were not associated with resistance to remdesivir. Functional characterization showed that mutations in the F gene led to a hyperfusogenic phenotype predominantly mediated by N465I. Additionally, recombinant wild-type-based MeV with the SSPE-F gene or the F gene with the N465I mutation was no longer lymphotropic but instead efficiently disseminated in neural cultures. Altogether, this case encourages further investigation of remdesivir as a potential treatment of SSPE and highlights the necessity to functionally understand SSPE-causing MeV.IMPORTANCEMeasles virus (MeV) causes acute, systemic disease and remains an important cause of morbidity and mortality in humans. Despite the lack of known entry receptors in the brain, MeV can persistently infect the brain causing the rare but fatal neurological disorder subacute sclerosing panencephalitis (SSPE). SSPE-causing MeVs are characterized by a hypermutated genome and a hyperfusogenic F protein that facilitates the rapid spread of MeV throughout the brain. No treatment against SSPE is available, but the nucleoside analog remdesivir was recently demonstrated to be effective against MeV in vitro. We show that treatment of an SSPE patient with remdesivir led to transient clinical improvement and did not induce viral escape mutants, encouraging the future use of remdesivir in SSPE patients. Functional characterization of the viral proteins sheds light on the shared properties of SSPE-causing MeVs and further contributes to understanding how those viruses cause disease.

Revisiting host identification of bacteriophage Enc34: from biochemical to molecular.

In our 2012 genome announcement (J Virol 86:11403-11404, 2012, https://doi.org/10.1128/JVI.01954-12), we initially identified the host bacterium of bacteriophage Enc34 as Enterobacter cancerogenus using biochemical tests. However, later in-house DNA sequencing revealed that the true host is a strain of Hafnia alvei. Capitalizing on our new DNA-sequencing capabilities, we also refined the genomic termini of Enc34, confirming a 60,496-bp genome with 12-nucleotide 5' cohesive ends. IMPORTANCE: Our correction reflects the evolving landscape of bacterial identification, where molecular methods have supplanted traditional biochemical tests. This case underscores the significance of revisiting past identifications, as seemingly known bacterial strains may yield unexpected discoveries, necessitating essential updates to the scientific record. Despite the host identity correction, our genome announcement retains importance as the first complete genome sequence of a Hafnia alvei bacteriophage.

Oxygen-dependent histone lysine demethylase 4 restricts hepatitis B virus replication.

Mammalian cells have evolved strategies to regulate gene expression when oxygen is limited. Hypoxia-inducible factors (HIF) are the major transcriptional regulators of host gene expression. We previously reported that HIFs bind and activate hepatitis B virus (HBV) DNA transcription under low oxygen conditions; however, the global cellular response to low oxygen is mediated by a family of oxygenases that work in concert with HIFs. Recent studies have identified a role for chromatin modifiers in sensing cellular oxygen and orchestrating transcriptional responses, but their role in the HBV life cycle is as yet undefined. We demonstrated that histone lysine demethylase 4 (KDM4) can restrict HBV, and pharmacological or oxygen-mediated inhibition of the demethylase increases viral RNAs derived from both episomal and integrated copies of the viral genome. Sequencing studies demonstrated that KDM4 is a major regulator of the hepatic transcriptome, which defines hepatocellular permissivity to HBV infection. We propose a model where HBV exploits cellular oxygen sensors to replicate and persist in the liver. Understanding oxygen-dependent pathways that regulate HBV infection will facilitate the development of physiologically relevant cell-based models that support efficient HBV replication.

Dinucleotide biases in the genomes of prokaryotic and eukaryotic dsDNA viruses and their hosts.

The genomes of cellular organisms display CpG and TpA dinucleotide composition biases. Such biases have been poorly investigated in dsDNA viruses. Here, we show that in dsDNA virus, bacterial, and eukaryotic genomes, the representation of TpA and CpG dinucleotides is strongly dependent on genomic G + C content. Thus, the classical observed/expected ratios do not fully capture dinucleotide biases across genomes. Because a larger portion of the variance in TpA frequency was explained by G + C content, we explored which additional factors drive the distribution of CpG dinucleotides. Using the residuals of the linear regressions as a measure of dinucleotide abundance and ancestral state reconstruction across eukaryotic and prokaryotic virus trees, we identified an important role for phylogeny in driving CpG representation. Nonetheless, phylogenetic ANOVA analyses showed that few host associations also account for significant variations. Among eukaryotic viruses, most significant differences were observed between arthropod-infecting viruses and viruses that infect vertebrates or unicellular organisms. However, an effect of viral DNA methylation status (either driven by the host or by viral-encoded methyltransferases) is also likely. Among prokaryotic viruses, cyanobacteria-infecting phages resulted to be significantly CpG-depleted, whereas phages that infect bacteria in the genera Burkolderia and Staphylococcus were CpG-rich. Comparison with bacterial genomes indicated that this effect is largely driven by the general tendency for phages to resemble the host's genomic CpG content. Notably, such tendency is stronger for temperate than for lytic phages. Our data shed light into the processes that shape virus genome composition and inform manipulation strategies for biotechnological applications.

Efficient Recovery of Complete Gut Viral Genomes by Combined Short- and Long-Read Sequencing.

Current metagenome assembled human gut phage catalogs contained mostly fragmented genomes. Here, comprehensive gut virome detection procedure is developed involving virus-like particle (VLP) enrichment from ≈500 g feces and combined sequencing of short- and long-read. Applied to 135 samples, a Chinese Gut Virome Catalog (CHGV) is assembled consisting of 21,499 non-redundant viral operational taxonomic units (vOTUs) that are significantly longer than those obtained by short-read sequencing and contained ≈35% (7675) complete genomes, which is ≈nine times more than those in the Gut Virome Database (GVD, ≈4%, 1,443). Interestingly, the majority (≈60%, 13,356) of the CHGV vOTUs are obtained by either long-read or hybrid assemblies, with little overlap with those assembled from only the short-read data. With this dataset, vast diversity of the gut virome is elucidated, including the identification of 32% (6,962) novel vOTUs compare to public gut virome databases, dozens of phages that are more prevalent than the crAssphages and/or Gubaphages, and several viral clades that are more diverse than the two. Finally, the functional capacities are also characterized of the CHGV encoded proteins and constructed a viral-host interaction network to facilitate future research and applications.

Genetic diversity of RNA viruses infecting invertebrate pests of rice.

Invertebrate species are a natural reservoir of viral genetic diversity, and invertebrate pests are widely distributed in crop fields. However, information on viruses infecting invertebrate pests of crops is limited. In this report, we describe the deep metatranscriptomic sequencing of 88 invertebrate samples covering all major invertebrate pests in rice fields. We identified 296 new RNA viruses and 13 known RNA viruses. These viruses clustered within 31 families, with many highly divergent viruses constituting potentially new families and genera. Of the identified viruses, 13 RNA viruses clustered within the Fiersviridae family of bacteriophages, and 48 RNA viruses clustered within families and genera of mycoviruses. We detected known rice viruses in novel invertebrate hosts at high abundances. Furthermore, some novel RNA viruses have genome structures closely matching to known plant viruses and clustered within genera of several plant virus species. Forty-five potential insect pathogenic RNA viruses were detected in invertebrate species. Our analysis revealed that host taxonomy plays a major role and geographical location plays an important role in structuring viral diversity. Cross-species transmission of RNA viruses was detected between invertebrate hosts. Newly identified viral genomes showed extensive variation for invertebrate viral families or genera. Together, the large-scale metatranscriptomic analysis greatly expands our understanding of RNA viruses in rice invertebrate species, the results provide valuable information for developing efficient strategies to manage insect pests and virus-mediated crop diseases.